N95- 19006 



1994 

NASA/ASEE SUMMER FACULTY FELLOWSHIP PROGRAM 



MARSHALL SPACE FLIGHT CENTER 
THE UNIVERSITY OF ALABAMA 


SCATTERING OF POINT SOURCE ILLUMINATION 
BY AN ARBITRARY CONFIGURATION 


Prepared by: Richard Solakiewicz 

Academic Rank: Associate Professor 


Institution and 

Department: Chicago State University 

Department of Mathematics 
and Computer Science 


NASA/MSFC 

Laboratory: 

Division: 

Branch: 


Space Science 

Earth Sciences/Applications 

Remote Sensing 


MSFC Colleagues: 


William Koshak 
Richard Blakeslee 
Hugh Christian 


xxxix 




Monte Carlo Simulations [1], equivalent medium approaches [2, 3], and methods based on 
Boltzmann transport theoiy [4, 5] have been applied to study the scattering of light produced by 
lightning in a thundercloud. The Monte Carlo approach is the most straightforward but other 
approaches may yield better insight into the phenomena. 

The cloud is assumed to consist of a uniform, homogeneous, random distribution of spherical 
water drops with an average radius of 10 pm and number density p = 100 cm -3 . In particular, 
light in the near infrared, X = 0.7114 pm, is of interest. The equivalent medium approach, based 
on methods due to Twersky [6, 7], yields the coherently scattered fields. For the ratio of the 
average separation distance between water drops to the wavelength, the scattering is considered 
to be almost totally incoherent [8]. During the time interval of a typical lightning event, the 
movement of the water drops is negligible. For fixed configurations of scatterers, the distinction 
between coherent and incoherent scattering is not very clear in the literature. The best explanation 
in the author's estimation is given by Foldy [9], who admits that for a fixed configuration of 
scatterers, all of the scattering is strictly coherent. A rather artificial distinction is made. The 
following definition is unambiguous: the square magnitude of an average field |<\|/>| 2 , \\f = E or H, 
is proportional to the coherent intensity; V = <|vj/| 2 ) - |<\|/)p i s proportional to the incoherent 
intensity. The function Fis approximately given in terms of the square magnitude of the radiated 
field from a scatterer at b, averaged over all possible configurations of the remaining scatterers, 

0) 

where dr s is a volume element in local coordinates r = r-b 

s s’ 

The problem of electromagnetic scattering of an incident plane wave by an arbitrary 
configuration of obstacles was solved by Twersky [10]. In this report, the results are extended to 
point source incidence corresponding to a Hertz dipole. Knowledge of the response of a fixed 
configuration of scatterers excited by a point source may provide insight to improve the accuracy 
of the values of bulk parameters for clouds which have been found using plane wave excitation. 

As in [3], we transform to the frequency domain; time domain solutions are recovered by a 
Fourier integral. A dyadic formalism is used throughout. We will employ the following notation: 
<|), w, and y will denote the incident, scattered, and total fields respectively. Lower case letters will 
correspond to single scattering (an object in isolation); upper case letters will be used for multiply 
scattered fields. Functional dependencies in brackets imply plane wave incidence; the use of 
parentheses is reserved for fields corresponding to a point source. Modification from plane wave 
to point source excitation is accomplished by operating on the plane wave results with a 
Sommerfeld-type contour integral representation [1 1] of a spherical Hankel function. 

^ Generally^ speaking, if a transverse wave has a direction f, its dyadic form will look like 
I - rr , where I is the identity dyad. Incident fields will have the form 


$[r, k] = (i - kk)(j)[r, k} 4>[r. k] = e**, $(r, r') = (i + ^)<|>(r, r<), <«r, r') = /£>(*|r - rfy (2) 


where the first argument refers to the point of observation and the second is the direction of 
incidence or the source location. The plane wave can be taken as the incident electric or magnetic 
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field. The function $(r,r') specifies the fields due a dipole located at r'. For an electric or 
magnetic dipole with directions p,p m , we have [3] 


E‘ = 



k 3 

4tccos 0 



H' 


k 3 


V/if-rfxp, 


r 'j)xp m H' m = -^r^(r,r')-p 


4k 


(3) 


where s 0 and |i 0 are the permittivity and permeability of the medium external to the scatterers. 
Here, E' and H' are the electric and magnetic field vectors due to an electric dipole; the subscripts 
m denote the corresponding quantities due to a magnetic dipole. We will solve for electric fields 
when given an electric dipole for the source and work with magnetic fields for a magnetic dipole. 
The remaining fields may be found using Maxwell’s eqs. 

A dyadic version of the Helmholtz surface integral representation is given by [10] 

“ = 4I7J dS(rf X $(r,r)] r • [V X v] - [v X $(r,r)f • [n x y]j = j$(r,r), v}. ( 4 ) 

where V operates on the variables associated with the vector r to a point on the surface of 
integration, T denotes the transpose, and n is the outward normal. In (4), \j/ = vj/[r,k] or \j/(r,r') 

A 

and u can be w[r,k] or u{ r,r') depending on the initial excitation. 


As r ~ qo, we can write 


4>(r, r) ~ (i — rr ^^{kr)e~ MT 

(5) 

in (4) to obtain 


i 

i 

it 

i 

II 

g 

Sr 

l 

(6) 


Single scattering amplitudes for both kinds of excitations may be defined in this way. 
Sommerfeld’s integral representation for /? 0 (1) is given by 


(l +‘|r)^ 1) ( A l r “ r D = ^{^f c )(T - (7) 

where r c ~ r c (0 c ,(p c ), 0 < (p c < 27 c, and 0 C starts at zero and goes to 0, - zoo, where 0 ; is in an 
interval which guarantees the convergence of the integral. Substituting in (4), reversing the order 
of integrations, and recognizing g from (6), yields the spectral representations 



^J^Q(f c )e^ r g[f c ,k], w(r,r') 


~ J dQ( r c ) e iktc ' r g(r c , r '). 


( 8 ) 


These representations are valid for r greater than the scatterer’s projection on r . 

Eq. (7) shows that a point source can be written as a superposition of plane waves. If the 
response of an object for plane wave excitation is known, its response to a point source can be 
obtained by superposition. Allowing the direction of incidence of a plane wave k to be complex 
and regarding exp[-/k -r'] in (7) as a phase factor, we find the relation 
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e~ ikT 'g[r , k| 



An integral operator like (9) will be used to modify plane wave forms to point source excitation. 

For spherical scatterers, explicit forms for the scattering amplitudes may be written in terms 
of Hansen's functions. The excitations may be written as 


}[ r > £ ] = Z N»(” + W^MUk, r)C_„(k) - r)B_(k)]. 

$( r > r ') = Z r ) M -™( k , r') + N 3 „„(k, r)Ni ra (*,r')], 


r>r': 


y==yy n*- ( 2n + 1 )(-i) m 

^ n(n + 1) 5 


m=1 m=—n 


( 10 ) 


where C and B are vector spherical harmonics. We use the definitions of these functions as given 
in [12] except that complex exponentials are used instead of even (cosine) and odd (sine) forms. 
The normalization of the spherical harmonics Y n m used here follows [11]. When r <r\ the 
corresponding vectors in the second eq. in (10) are interchanged. 

The interior and radiated fields have the forms 


( 11 ) 


V»[ r > k] = Z r)a”[k] + NLlX, r)b"[k]J, 

«[ r > k] = X «0c[k] + N 3 „<X, r)d”[k]|, 

V, ir,r') = £ {MU*,r>:(r') + NL„(AT, r)b"(r')}. 

"( r * r ') = Z {ML(^, r)C(r') + NUK, r)d:(r')}, 

where k = £(e 0 ,p 0 ) and K- £{e,p) is the interior wave number. In general, the transition 
conditions at an interface, where the normal is denoted by n, require that 

n x (<j> + w) = n x vj/-i, n x V x (<j> + */) = Bn x V x \j/. w , (12) 

where B = p 0 /p if $ represents an electric field and e 0 /e if <ji is taken to be a magnetic field. 

For a sphere of radius a , the scattering coefficients are given by 

a„”[k] = rf&WZC. 4k)a„, a”(r') = QTMl_(A,rX, 

b:[k]=/”- i x»+i)e;B_„(k>„ J K(r’) =Q:^_jky)b n , 

c:[k]=/"V»(»+i)e;c m „(k) c „ > ^ 

<i:[k] = /"-■>(« +i)a m B_„„(ky„, n:(r') = e; o.rx. 

where the constants a n , b n , c n , and d n are the appropriate set of Mie coefficients found in [3], 
Series for the scattering amplitudes may be obtained by substituting large kr forms for the 
Hansen's functions in ( 1 1 ), 


(13) 


g[r,k] = ^ + l)&»Jc'„C m „(r)C „„(k) + <B„„(f)B_„, ,(£)], 

f (r, r-) = £ V«(« + V'-"Q:[c„C m Ji)M'_Jk, r-) + r')]. 


(14) 
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and may be used for the electric or magnetic fields given the same was chosen as the incidence 
after multiplication by appropriate constants. The remaining fields may be found by first 
interchanging those M's and N's in u or C's and B's in g which are functions of r. Solutions 
u and g for magnetic fields w hen $ is taken to be an electric field are then obtained by 
multiplication by y = -i^e 0 /[i 0 ; use 1/y for the remaining situation. 

Consider a fixed configuration of N scatterers. Each has its "center" located at 
, s = 1, 2, . . . , N. The total field for plane wave excitation may be written as 

$[r, k] = $[r, k] + U[ r,k] fi[r, k] = £ u{ r„ kje* 1 ’. ~ /<®(Ar)G[r, kj r, = r-b„ ( 1 5 ) 

where U s is the scattered field from the scatterer at b s as if it was located at the origin. Forms 
similar to those obtained for an isolated scatterer. 


k] = 2^ f *G s [r„ k] ~ k], G,[r,k] = |(i -ff^-^.t/Jr.kjl ( 16 ) 
where r ~ r as r ~ co, are also valid here. Taking phase differences into account, we have 


(17) 


G[f,k] = 2^" fK ^[f,kl 

s 

Twersky [10] obtained a coupled set of integral eqs. for GJr,k] in terms g, 

$[?.£] = &[f, k] + ^ £ J *J(i- ) r c ] . G,[r„ k} b„ =b,-b„ (18) 


Recall that G s is related to U s , which is written in local coordinates r s . Operating on (18) with 
the integral JQ^k^ex p[-/k-(r'-b,)] yields 


G,(r, b;) = &(r,b() ?„]• G,(r c , b;), b; = r'-b,. ( 19 ) 


When the particle separations are large compared to wavelength, the integral in (19) may be 
evaluated asymptotically. A convenient method [10] to do this exists whenever the kernel, aside 
from the exponential, can be expanded in a series of vector spherical harmonics. Iterating while 
retaining an appropriate number of terms in the asymptotic series will generate a consistent 
asymptotic expansion for G s . Exactly the same expansion as given in [10] will work for (19) 
except that G[] is replaced by G(). Retaining only the leading term in the integral, gives 

G(r , b() ~ £(?, b;) + £ b J • G,(b„ , b;). (20) 

A set of algebraic eqs. may be obtained by substituting (14) and 

G,(r,b;) = X [c„(f)C(„,(b;) + B„,„(f)D* n (b')] (21) 

into (19) and using orthogonality. Similar to [10], we obtain 
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( 22 ) 


CL = Mn+T^gSSLjk, bf ) + rtn + 1)Qk£ Z [^(b «)Q* P + Fi(b,)D^l 

. *** Q’P 

^ + ty- n d s n QZW_ mn (k, b') + */i + 1)2K£ Z [-^(b,)a + E^b^l 

q,p * 

£ -"(M = ^ JdO(4)«^ k 'C-_fc) • cjr,) = -L J dQ(r c )e^B_Ji) ■ B„(r 0 ), 
F -”"(b«) = ^ J dn(n,)«"'' b "C4r e ). B„(f c ) = -ij rfn(f t ) i ;*--‘.B_(f t ) • Cjf,), 


where the scalars fand Fean be evaluated in terms of special functions. 

In a cubic cloud 10 km on a side, there are approximately 10 20 drops. For such large systems, 
a simple approach as suggested in (20) may be appropriate. Sufficiently far from the source, each 
drop in a small sub-volume of the cloud experiences essentially the same initial excitation. All of 
the drops in this sub-volume can be considered as a "compound” scatterer. Its scattering 
amplitude may be determined by (17). It will then be possible to consider fewer particles making 
up the cloudy Presumably, the procedure of obtaining G from G, considering G as a 
"compound" G, and working with a distribution of "compound" scatterers can be repeated to 
obtain a result which is amenable to numerical computations. 

The author would like to express his appreciation to W. Koshak, R. Blakeslee, H. Christian, 
and D. Phanord for their help and acknowledge the assistance of F. Six, L. Freeman, G. Karr, S. 
Novy and T. Shurtz. 
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